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1. Introduction 

Seemingly simple quantum mechanical problems like the traversal of potential barriers 
by wave packets are surprisingly rich of mathematical, technical and physical details. 
The 1-dimensional case is interesting in its own: electron transport through barrier 
junctions, the physics of light and wave guides, particularly of the transmission through 
Photonic Band Gaps, the time of arrival of wave packets and the ensuing paradoxes of 
non-locality and super-luminal tunnelling are examples. 

Because of its strong overlapping with the standard theory of 3-D scattering, usually 
the s-wave radial problem for the spherical shell potential is studied in the literature 
as an example of the role of the potential barriers in connection with the evolution 
of decaying states. These states obey specific boundary conditions (BCs), namely the 
vanishing of the wave function at the origin and a purely outgoing BC outside the 
barrier. However the corresponding BCs adopted in the 1-D scattering make it different 
also from a mathematical point of view, providing a different setting for the related 
resonant (Gamow) states. 

On the other hand, practical calculations face important limitations and the effort 
of pushing the analytical calculation as far as possible is always rewarding. Significant 
progress has been accomplished [1] for a Gaussian packet impinging on a square barrier, 
where the Gaussian structure was exploited, and [2] for the shutter problem, where the 
contribution of the resonant states in a double square barrier was worked out. Relevant 
resonances are mostly expected in systems of two or more barriers separated by a gap. 
Actually they generally occur in any simple barrier, and the plain square barrier is most 
tractable, fully representative for many theoretical purposes and devoid of bound states 
or other unessential features. 

In this paper we consider the time evolution, in the presence of a square barrier, 
of an initial plane wave with support on the whole 1-D space, aiming to work out the 
contribution of the resonant states. These states have a central role, as long as they 
correspond to complex momenta (or energies) which are poles of the S-matrix, and 
dictate the most convenient type of BCs to be used in the problem. 

Our calculation extends the results of shutter problem, where the initial plane wave 
occupied the space axis at one side of the barrier and was specially suited to study details 
of the propagation of the wave front. We provide a comprehensive result containing also 
the evolution of the segments of the initial plane wave lying inside and at the other side 
of the barrier. 

Our main goal however, is to pave the way for further uses since superpositions of 
infinite plane waves build up any desired initial wave packet. In particular, the compact 
result obtained for the Gaussian wave packet pj could be re-derived pinpointing the 
contribution of each single Fourier component. 

A byproduct of our calculation is to provide an explicit example where we can study 
the properties of completeness of the resonant states in different regions of the square 
potential, an issue properly dealt with [3J in the frame of Rigged Hilbert Spaces, which 
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will be addressed elsewhere. 

In this paper we follow the approach used by Peierls and Garcia- Calderon [1], 
namely Laplace-transforming the time evolution equation into a second order linear 
differential equation (SOLDE) in one variable, expressing its solutions in terms of the 
Green function (GF) with resonant boundary conditions (RBCs) and then undoing the 
Laplace transform. 

The first step is done in Section 2, where a brief discussion is made of the use of the 
Green's method for the solution of SOLDEs when the BCs required for the GF and for 
the solution are different. This point is relevant for the issue of the completeness of the 
Gamow states. The GF with RBCs features a simple structure of isolated (resonance) 
poles, in terms of which the (p-dependent) Laplace-transformed wave function is worked 
out in Section 3 previously to the second step. Then the explicit calculation of the 
inverse Laplace transform leading to the sought after t-dependent solution is carried out 
in Section 4. 

As a check of this cumbersome analytical calculation, the t = (Section 5 of the 
paper) and the t = oo (Section 6) limits of the time-dependent solution ip{x,t) are 
calculated: at t = one must recover the initial plane wave and at t = oo a suitable 
stationary solution must be reached. Then the conclusions are drawn in Section 7. 

Some notations and a number of calculations and technical details are deferred to 
the Appendices in an effort to make the paper more self-contained. 



2. The evolution equations 

We consider the 1-D time-dependent Schrodinger equation 

{ih^^-H)^(x,t) = (1) 

where the Hamiltonian H = / 2m + V{x) corresponds to the square barrier potential 
V{x) = 6{x)9{L — x)V and the solution satisfies the initial condition ip{x, 0) = ipo{x) = 
e'''^ , A; > . 

The Laplace transform ipls) = dt e~^^ilj{t) on the time variable, brings the 
parabolic partial derivative differential equation (1) to the simpler SOLDE 

' +p'~-^V{x)]i^{x,p)=iae"'^ (2) 



dx^ h 

where a = 2m/h and p'^ = ia s (Appendix A). The homogeneous equation L^ip = 0, 
namely the time-independent Schrodinger equation for the potential V{x), has a variety 
of solutions according to the BCs we choose at x = and x = L depending on the 
physical problem. Besides the usual scattering in and out solutions for continuous 
and real E = p^ / 2m > 0, one has the resonant solutions satisfying the homogeneous 
outgoing RBCs 

5x^U=o= -ipV^(O) , d^ijj\^=L = ipi'iL) (3) 
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(the ingoing reversed sign of p is not considered here for the reasons outhned later on). 
These solutions exist only for a dcnumcrable set of isolated values p„ of p lying in the 
lower half complex plane (Appendix A). 

A solution to the inhomogeneous equation (2) for e can be explicitly written 
down for each region of V{x) : 

Mx,p) ^Be-'P'^ + ^^e"'^ (x<0) 

iPn{x,p) =Me^^'" + iVe-'^''" + ^j^e'""- (0 < a; < L) (4) 



where the amplitudes A, B, M and are functions of p completely determined by the 
matching conditions at a; = and x = L and have a common denominator of the form 
D(p) = 0^e'^'^ — ©^e~'^''^ (notation in Appendix A). Out of all the terms of the general 
homogeneous solution, the choice of e'^^ (for a; > L) or e~'*'^ (for x < 0) is dictated by 
the behaviour of the solution for x — > ±oo , since a small positive imaginary part in p 
is supposed when performing the inverse Laplace transform leading to ip{x,t) , t > 
(Appendix D). 

However, this expression of the solution is not suited to perform the inverse Laplace 
transform back to ip{x,t) because of the nontrivial analytical form of the amplitudes 
above. Instead the GF approach lets us to express ilj{x,p) as a sum of isolated pole 
terms plus other simple terms easier to deal with. To this end, (2) together with the 
Green equation for G{x,x',p) can be written in the form 

L^^(x,p) -iae^'^^ =0 (5) 

L^G{x, x',p)- d{x -x')=Q , (6) 



where G{x,x',p) is required to obey RBCs, namely 

da;G{x,x',p%=o^ -ipG{0,x',p) , da;G{x,x',p%=L^ ipG{L,x',p) (7) 

It is expected that G{x, x',p) = will have poles at p = p„ , where the homogeneous 
equation L^'ip — with RBCs has non-trivial solutions. 

The equations (5) and (6) may conveniently be given the short-hand notation 
^ — and r = respectively. Then the Green method uses the integral equation 

Jq dx [ijjr — G'^] = to obtain 'ijj{x,p) in terms of G{x,x' ,p). Only if t[j and G obey the 
same homogeneous BCs, the surface terms in this integral cancel out and one obtains 
the familiar result ipn{x, p) = Jq dx' ia G{x' , x, p) e^^^' . 

However the BCs for ilj{x,p) (involving ip and d^ip both at a; = and at a; = L) 
stemming from (4) are different from (3) and non-homogeneous. In that case one obtains 
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a modified expression for tlie Region II, namely 

il)ii{x,p) = / dx' ia G(x', e^'^^ 
Jo 



G{L,x,p)e'''^ -G{0,x,p) (8) 

p + K p — K 

For tlie external regions / and ///, tlie matcfiing of at x = and x = L determine 
tlie coefficients B and A respectively, so that 

i>iix,p) = iae-'P'' dx'iaG{x',0,p)e'^''' 

Jo 

" 0, j9)e''=^e-'P" - -^GiO, 0, p)e-'^'" 



p + k p — k 



and 

i,ni{x,p) = iae'f("-^) dx' iaG{x', L,p) e""""' 

Jo 



+ I^e'- (9) 



p + k p — k 

1'^ ^ip{x-L) ifcL _|_ ifcx 



For comparison, only the last term in the r.h.s. of (8), the last three ones in (9) and the 
third one in (10) appear in the shutter problem. 

The crucial advantage of the method is that almost all the Green functions involved 
in the equations above can be expanded as a sum of terms which are simple poles in 



3. Analytical structure of the p-dependent solution 

General theorems ^ and the explicit analytical derivation (Appendix B) of G{x,x',p) 
show that \G{x,x' ,p)\ ~ l/\p\ — ^ as |p| — oo in the complex plane for almost any 
X G [0, L] and x' G [0, L] . Then the Mittag-Leffler theorem tells that 

G(x,x',p) = E^^^^^ (11) 

n P~Pn 

The exception is for G{0,0,p) and G{L,L,p), the modulus of which grows as \p\ for 
IpI ^ oo in the lower half complex plane p and requires some special care (Appendix 
C). 

The residues in the r.h.s of (11) can be easily computed |1] and one finds 
Gn{x,x') = Un{x)un{x') / Nn , whcrc the functions m„(x) belong to the denumerable 
set of the resonant solutions satisfying L^. m„(x) = with RBCs 
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and arc suitable normalization factors (Appendix A). 
The inverse Laplace transform 

I /-c+ioo I r+oo . p2 ^ _ 

Hx,t) = — / dse"V(a^,p(s)) = 7^ / dppe '^m^ipix^p) , (12) 

ZTTl Jc-ioo lum J -00 

written as an integral over the real momentum variable p (Appendix D), leads one to 
consider the pole expansion of p ilj{x,p) in the integrand, which can be worked out 
for each sector /, // and ///. We report here only the result (Appendix E) for the 
sectors // and /// (sector I is similar to sector III because of the symmetry of both 
the potential and the initial condition): 

p Un{x) rL 



pi^n{x,p) ^ia^^—'^f dx' u^{x') e^'^' 



p-t-h, n P Pni^^Pn -'*n 

- a--^GiO, x,k)+a^ ,-,„k-p. ^''^ 



pM^.p) ^i«e'^(--^)^ ^ f^^dx'unix')^^^' 

n P-Pn Nn Jo 

2 

- a^GiL, L, 0)e^P(^-^)e''=^ + -^jG{L, L, -kyp^^-^^)^^^ 

aY P^ ^ ^ ^l(^) ^ip(x-L)^ikL 
n P-PnPnk + Pn 

+ ^ Pn 1 i^n(OK(^) ^ip(.-L) ^ ^_^G(0, L, A;)e'^(--^) 

n P~Pnk—Pn k—p 

i« ^^)(e''=^ - e^^'^^-^^e^^^) (14) 

One recognizes the terms of the shutter problem in the third row of (13) and in the 
fourth row of (14). 

4. The time-dependent solution 

For each of the terms above, the integrals stemming from (12) can be brought to the 
form of an integral representation of the (complementary) error function erfc{z) so that 
their inverse Laplace transform can be carried out thoroughly (Appendix F). We obtain: 

Unix) 



+ ikG{L, X, -^)e^^^ [B^^] - ie'^^ ^ ^n{L)u^^x) ^^..^ 

n 1^ Pn -'^n 
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- ikGiO, X, k) [S^^] + i E ^-(QK(^) ^gn^ (15) 

n Pn -''n 

where the factors in square brackets embody the time (j = t / 2m) dependence of the 
solution and are the result of the integrations over the momentum p. Explicitly 



— 


— er/c(i-\/iT) 


B'' = 


— e^ 


erfc{-iViTk) 


Bl! = 


— e~ 


-^-rpi erfc{iViTpn) 


S'' = 


— e" 


-'^^^ erfciiVwk) 


all _ 


— e" 


-i-rpn er/c(i\/irpn) 



Likewise 



n JO 



ikL ikL/ 

- iG{L, L, of— [B',''] + iG{L, L, -kf— [BLi'] 



(16) 



n Pn rv ^ Pn 

- ikGiO, L, k) [S^^^] + i E MO)ML) ^^ni^ 

n Pn -''n 



where 



e' 4r (p^e^" er/c(y„) + —=) 

J-KT 



Xlll 



r777 e'4 X-L .(x^ 

-On = : e -It 

Sii^ = - e^'-^ik^yl erfc{y^,) + -^{-k + ^)) 

y/TTT It 

■(x-L)'^ n 2 e'4 a; — L 

Bi'' = - e^^ipy'^ erfciy^) + ^(pn + -—-)) 

yTTT It 

S^ii = - e^^^e^'' erfcivk) 

Si!' = - e'^e^" er/c(i/„) 
^ = - e-'^'^' erfc{-iV^k) 
Ci^^ = - e"'^*^' erfciiV^^k) 

Ci'' = - e'Te^^-. er/c(y_,) 
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Cf ' = - e'^eyl erfcivk) (18) 

We have introduced the variables yq = c^'t (4r)^5 — L — 2rg) for q = pn, k, —k. In 
both equations the terms [S] are the ones arising in the shutter problem. The particular 
values of the Green function involved in eqs.(15) and (17) are calculated in Appendix 
B. 



5. The short time Umit 



The t ^ Q limit is interesting both as a check of the calculation above and for the 
study of the scattered wave at short times. Here we aim only to recover the initial wave 
function il!{x,0) — e^''^ at t — 0, which must happen in each of the sectors /, // and 
III. The sector I similar to so we refrain from calculating it. 

Careful use of the limiting values of erfc{z) and/or w{z) — e'~^^er/c(— i^;) for 2; — > 
and for 2; — cxo in different directions of the complex z plane must be made. Notice 
that Uq tends to 00 in different directions for different q. Also the properties of the set 
of resonant functions m„(x) as a basis of the space of solutions are crucial (Appendix 
G). 

5.1. Region II 

For r ^ the factors [B^% [B^J], [S^^] and [S^^] tend to the value -1, so that the [B] 
terms in (15) yield 

Vn 

{L)Un{x) 



ikG{L, X, -ky^^ - le'^^ Yl 



k + Pn Nn 



— \ke^^^ 1 Un{L)Un{x) ^ .^i^j^, Pn Un{L)Un{x) 

^-k-pn Nn nk+Pn 

^_^^.kL^ ML)Un{x) (19) 



and hkewise for the [5"] terms: 

Pn Un{0)Un{x) 



ikG{Q,x,k) -iY^ 



k-Pn Nn 

^k~Pn Nn '^k-pn Nn 

^.^MOM^^O (20) 



The limit of the factor [A^^] is not readable in (16) since it seemingly blows up for 
r — > 0. Actually, the changes of variables done to perform the integration leading to 
this form of [A^J] become singular in this hmit. The right result can be reached by the 
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following reasoning: For r — > 0, an increasingly large range of values of p (as compared 
to Repn) does contribute to the integral dp e~^'^^^ p / {p — Pn) ■ Therefore its limiting 
value (tt / ir)^ does not depend on pn and becomes a common factor outside the first 
sum in (15). Then one recovers the sum J2un{x)un{x')/Nn = S{x — x'), as expected 
from the t — > limit of the Green function G{x, x', t). 
Therefore, for < a; < L , 

limt-^o i^ii{x, t) = e'*^"^ = il^oix) 

as required. 

Asymptotic expressions for erfc{z) and w{z) can be used to obtain approximations 
to the form of '^ii{x, t) for small values of t, a task that will be faced elsewhere. 

5.2. Region III 

The limit must be studied term by term. The integral dp e'^'^^^&^^'^'^^p/ [p — Pn) 
has now an extra exponential factor which has a regularizing effect. As a result, the 
limit can be taken directly in [Al^^\ as given in (18): the term e^^erfcd/n) vanishes and 
(7rT)~5e'' 4t^ approaches a distribution concentrated in x — L, namely 6{x — L). Thus 
it vanishes for x > L. 

An analogous reasoning shows that [-Bq^^], [B[^^] and [B^J^] yield derivatives of 
5{x — L) and also vanish for x > L. 

The factors [S^^^], [S^^^], [C|^^] and [Ci^^] vanish exactly, whereas [Cf^^] and [C|^^] 
tend to —1. From (17) we then see that the final result is that also 

limt-^o ipiii{x, t) = e'''^ = ipoix) 

for X > L, as required. 

The solution ijji{x,t) at the left of the barrier will have a structure similar to (17) 
and the same limit above. 

6. The large time limit 

For r — > oo we see that also — > oo as before, but in still different directions of the 
complex plane for the different q. 

In the internal Region II, the factors [A^^], [B^^], [B^^] and [5^^] vanish, whereas 
[S^^] -2e-'^^\ Therefore, for < x < L, 

limt^oo i'liix, t) = 2ik G{0, x, k)e~''^'^^ 

= 0;"(a;)e-'-^'=* (21) 

that is the stationary scattering in solution, where = k"^ / 2m (see Appendices B and 
H). Thus the infinite plane wave evolves into the same final state of the shutter initial 
condition. 
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For the external Region III, the only non- vanishing factors are 

(22) 

Then the terms corresponding to C|^^ and C|^^ in (17) cancel each other and, again, 
one is left with the same asymptotic solution of the shutter 

Umt^^ ^in{x,t) = 2ik ^(0, L, A;)e''=(^-^)e-'"^' 

= 0t"(a;)e-^^'=* (23) 

The finite time behaviour is thus made up of transient modes which quickly 
disappear leaving (for k > 0) the scattering asymptotic solution with outgoing BC 

clIj X —— L. 

7. Conclusions 

A solution for the time evolution of an infinite plane wave in the presence of a simple 
square barrier has been worked out for each of the regions of the potential and Section 
4 is the main result of this paper. Among other terms, this solution contains a sum of 
explicit analytical contributions corresponding to each of the (infinitely many) resonance 
poles. As in previous works in related problems [1] [2], only the location of these poles 
needs to be obtained by numerical methods. 

In the solution obtained the shutter terms have been pinpointed. Similarly, the 
contributions to the evolved wave function coming from the segments of the initial wave 
function lying inside and at the right of the barrier can be identified. 

As the main application of this knowledge we envisage the possibility of studying the 
enhancement or the suppression of the transmission of the single Fourier components 
of any realistic wave packet. This should provide new detailed insight on interesting 
phenomena like the super- luminal tunnelling |6j, the breakdown of energy conservation 
[7j by transient interference in wave packet collisions with barriers [8] or the rising of 
forerunners. 

Other side results of this work are more mathematical and preliminary for further 
study. The exact analytic expression for the GF of the square barrier subject to RBCs 
has been worked out. Some particular values of it are involved in the time-dependent 
solution obtained above, but its knowledge is also useful for computing in each region 
the extra terms in general resonance pole expansions of wave functions in order to 
implement, with an explicit example, the related question of the completeness of the 
set of the resonant solutions We have relied on the adoption of different boundary 
conditions for the GF and for the solutions of the (Laplace transformed) evolution 
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equations. The link between this point and the properties of completeness of the 
solutions will be addressed elsewhere. 

On more physical grounds, the explicit solution obtained for finite time is useful for 
deriving approximations valid for short times, hence for the study of transient structures 
and forerunners. An immediate result of the work is that the large time limit yields the 
same stationary solution of the shutter. This again shows that the resonances contribute 
only to transient structures of the scattered wave. 
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Appendix A. Resonant solutions 

We adopt units such that h — 1 and define the differential operator 

L, = [—+p^-2mV{x)]. 

A piece-wise general solution for the homogeneous SOLDE u{x) = is 

u{x,p) = e{-x)Be-'P'' + e{x)e{L - x)(Me'f'^ + Ne-'P''') + e{x - L) Ae^^^ , 

which has already built-in the RBCs (3). Choosing B{p) arbitrary, the matching 
conditions at x — and at x — L yield the amplitudes 

M{p) = - ^B{p) 

m = - ^mp) 

A(p) = - ^(©e-^®^ - ee-'^^)B(p) , (A.l) 

where p' = — 2mV , ® = p + and Q = p — p' , plus an extra condition on p 
which is proportional to 

©2 g-ieL _ e2g-ieL _ (A.2) 

This transcendent equation has a denumcrablc set of solutions Pn lying in the lower 
complex plane p . The common denominator D(p) = 0^e'^''^ — ©^e"'^'"^ in the amplitudes 
of (4) is just proportional to the l.h.s. of (A.2). One can check that ii Pn is a solution, 
then — p* is too, so that these values are in symmetrical locations with respect to the 
imaginary axis. We let the label n take integer values {n ^ 0), with the growing 
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positive n indicating the Pn with growing real positive part and = — j9* . One finds 
that —J < argpn < and vr < aryp^n < x- 

The simplest writing of the solution (24) is obtained for B{p) = —2p', namely 

Unix) = e{-x){-2p'Je-'P"^ 

+ e{x)e{L - x)(e„e^P"" - ©^.e-'^-") 

+ e{x - L) (©„ e-^®"^ - ©„ e-^®"^) e^^'"^ , (A.3) 



Notice that the solution (4) should yield the outgoing resonant ones, which exist 
only for the momenta Pn , when we let the terms involving the initial condition e^''^ 
(namely the particular solution of the inhomogeneous equation) vanish. In fact, one 
sees that for p ^ Pn , these terms belong to the regular part of (4), with corresponding 
inhomogeneous BCs, while the pole part yields the resonant solution. 

The crucial role of the BCs in solving the Schrodinger equation [10] is behind the 
fact that different BCs like the RBCs, obeyed by the resonant solutions and by the Gf 
(causing its resonance poles), on one side, and the ones characterizing the scattering 
solutions (whose amplitudes display the S-matrix poles) on the other, lead to the same 
poles. 

The residues Cn{x,x') = Un{x)un{x') / Nn obtained in (11) correspond to the 
following choice of (complex) "norm" |1] 

Nn = KuliO) + uliL)) + 2pn dx ul{x) (A.4) 

JO 

which takes the value Nn = —8mV{pnL + 2i) for the solutions (A.3). This "norm" is 
somehow related to the usual one in Hilbert space: when applied to ordinary square- 
integrable wave functions ip^x) G C^, replacing the squares by the square modulus and 
letting and L recede respectively to —oo and to +oo, equation (A.4) approaches the 
usual norm (up to a factor). For the resonant solutions however, diverges in this 
limit . 

Appendix B. Analytical Green function 

The analytical solution to the Green equation LrcG{x,y,p) = 6{x — y) for x and y in 
the interval [0,L] , and obeying the RBCs (7), can be directly constructed: 

G{x,y,p) = 7^j^ { 2mV (e^p'(.L-{x+y)) ^ ^-ip'(L-(x+y))^ ^^ -^^ 

zp JJyp) 

- ®h-'P'^{e{y - a;)e'P'(^-^) + e{x - y)e'P'^^-y^) 

- ©'e'^'^(^(2/ - a;)e-'^'(^-") + e{x - y)e~'P'^''-y^) } , 



where the symmetry x ^ y is explicit and the zeroes of D{p) give the expected simple 
poles. 
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G(L,a;,p) = -^(Ge'^''" - ©e-'^'") 
D{p) 

2n' 

G(0,L,p) = -i ^ 



Dip) 



G(0,0,p) = G{L,L,p) = 7Tf^(©e'P'^ - ©e-^''^) 

D{p) 



-1 



G(x, y, 0) = , ; {Cosh[V2mV(L - (x + y))] 

+ 9{y - x)Cosh[V2mV{L - {y - x))] 

+ e{x - y)Cosh[V2mV{L - {x - y))]} 
1 



G(0, 0, 0) = G{L, L, 0) = Goth[LV2mV] 

\j2mV 

5,G(0,0,p)|,-o^ 1 3 + Co^^^2L^ 

The limit |p| — > c)0 can be directly read out in the above expressions. It can be 
seen that in all the cases the Green function vanishes as Xjp or faster in any direction 
of the complex plane, with the only exception of (^(0, 0,p) and GiL^ L,p) , which grow 
as p in the lower half plane, though they still decrease as 1/p in the real axis. 

The particular cases G'(0, x,p) and G{L, x,p) are related to the scattering solutions 
(see Appendix H). From the integral formula dx [(p {LxG—6{x—x'))—GLx4>] = 
and the use of the BCs for (p and the RBCs for G , one obtains fU] 



G(0,x,p) = — v^./^C(x) (0<x<L) (B.3) 
2p V m 

G(L,x,p) = ^v^,/^C(a:) (0<x<L) (B.4) 
2p \ m 



Appendix C. Pole expansion and substractions 

Whenever the Green function G{x,x',p) vanishes as 1/p in the limit \p\ — > oo, the 
Mittag-Leffler theorem applies and the validity of the pole expansion (11) is assured. 
This is not the case for G{0,0,p) and G{L,L,p) , which grow as p in the lower half 
plane, and a substraction technique must be used j5j. This is performed by using the 
contour integral 
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where the contour T = {Co, Cp , C„, Cg} consists of small closed paths counterclock-wise 
encircling the poles of the integrand (namely Co around the double pole ai z = 0, Cp 
around p and C„ around each pn) plus a large circle Cg clock- wise encircling all these 
poles. The factor z~'^ has been introduced to assure that the integrand goes as even 
for G{0, 0,p) and C{L, L,p) so that the integral over Cg vanishes when the radius of the 
circle grows to oo. 

From (C.l) one readily obtains (see also [T2] ) 

Gix, x\ p)=p'f: 1 ^M^KM + [G(x, x', 0)+p dpG{x, x', p) |,=o] (C.2) 

n NnPi P-Pn 

When the pole expansion (11) holds, (C.2) becomes a trivial identity, whereas in the 
case of G(0, 0,p) and G{L, L,p) a holomorphic (linearly growing with \p\ ) part remains. 



Appendix D. Inverse Laplace transform 

The inverse Laplace transform in the variable s in (12), involves an integration in the 
complex s-plane along a line parallel to the imaginary axis with Re c > 0. In the p-plane 
(recall p = ^/i^/2ms), this path translates to a hyperbola-like one with asymptotes in 

the positive real and imaginary axis. For t > , the factor e~^2^* assures that it can be 
deformed into an integration from +oo to — oo along (and slightly above) the real axis 
if the integrand has poles only in the lower half plane (and on the real axis). For t < , 
the path can be closed along a quarter circle of large radius in the 1st quadrant, thus 
enclosing a region without poles and giving ip{x,t) = 0, consistently with the Laplace 
method and causality [5]. 

Appendix E. Pole expansion of pijj{x,p) 

We rewrite (8) as 

M^,p) = ia ("^ dx' fi'{x\x,p)e''^' -ae'''^fi'{x,p)-afi\x,p) (E.l) 
Jo 

where 

fi\x\x,p) = G{x',x,p) , fi\x,p) = and M\x,p) = ^^^^ . 

p + k p — k 

According to the asymptotic behaviour of the fi\p) {i = 1,2,3) above we consider 
the contour integrals 

= — fdz f!\z) (E.2) 

ZTTl JFj Z — p 

where the values ni = 0, and n2 = = 1 are assigned so that the integrand of (E.2) 
behaves as z~^ for large z . The contours Fj include circles around the poles in the 
respective integrand (namely p , pn, k or —k ) plus the large circle Cg , as in Appendix 
C. For i = 1, equation (E.2) yields f(\p), which must be multiplied by p later on. 
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whereas for i = 2,3 one directly obtains p fl\p)- Collecting these results and using the 
pole expansion of the Green function, (13) is obtained. 

The derivation of (14) follows the same lines. Here we define 

^nu^\ ni^< T tlnl^^ G{L,L,p) ^jjj, ^ 0(0, L,p) 



fl^\p) = G{x',L,p) , fi^^p) = 11^^^ and f^'\p) 



p + k p — k 

and the powers required are ni = 0, n2 = — 1 and ns = 1 . Notice that the negative 
power = —1, needed by the asymptotic behaviour of G{L, L, z) , introduces an extra 
pole at z = and is related to the discussion of App.C. 

Appendix F. er/c f^zj-related integrals 

In the computation of ipn{x,t) and ipiii{x,t), the following types of integrals arise: 

J(x) = — dp p e-'"^''e'(^-^)'' (F.l) 
ivr J-oo 

1 r+oo 1 „ 

Ux, q) = — dp e-'^P qK^-l)p ^ _2 m(x - L, q, t) (F.2) 

ITT J-oo p — q 

h{x, q) = - dp e-'^v\^(^-L)v = -i|_/o(a;, q) (F.3) 

ITT J— oo p — q ox 

1 r+oo r)2 . 2 ■/ r^ d 



hix, q) = - dp^— e--P e^^^-^^P = (-i— )^/o(x, q) , (F.4) 
ITT J-oo P — q ox 



where 



/• + 00 

M{x,q,t) = — dp 

ZTTl J-oo 



_X r + OO I .42 



p-q 

is the Moshinsky function, often found in the literature. 

They can be related to the integral representation of the error function [13 



2 i r~r^ e 
w{z) = erfc{—iz) = / du Imz > (F-5) 

7r J-oo u — z 

so that 

M{x,q,t) = -e'-^t^ w{ie 'Y^(^-— 9))- 



When converting the exponent irp^ into v?, namely u = ^/irp , an integration along 
a line running between ±\/i oo results. In the case q = Pn, rotating this path by an angle 
7r/4 towards the real axis crosses the pole \/irp„, contributing a residue. For q = ±k, 
a small negative imaginary part is supposed (see Appendix D). The case Imz < can 
be tackled by using the property w{z*) = w*{—z). 

The results of the integrals above, as given in Section 4, are 
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= I{x) 






= /2(x, -k) 


Bl!^ 




= hix.Pn) 


giii 




= /o(x, fc) 


qIII 




= Io{x,Pn) 


B'' = 




= h{L,-k) 


S'' = 




= Io{L,k) 






= Io{x, -k) 






= Io{x,k) 


Bl! = 


oil 





2 

Appendix G. t — 0, oo limits of erfc{yq) 

— i ^ 

Being Uq = ^^^{x — L — 2rg) , where r = ^ , we see that, for t , 

_ii /^/ r\ -i- 

l/g~e iJ — {x-L)^e 4.00, 

regardless of q being equal to ±k, Pn or p_„ . 
For t —>■ oo one has 



t 



where f < 0„ < ^ , -f < 0_„ < , 0fc = ^ and 0.^ = -| . 
The following asymptotic formulas [13\ apply: 

71 

liniz^ao erfc{z) = \ci^) z\ < — 

2 11 °° C 

lim^^oo e"" erfc{z) ~ ^-(1 + ^) ^ l^-^d z\ < 

Only the case 4>k = ^ is out of their range, in which case the relationship 

qvI erfcijjk) = 2e^^ - e"^^ erfc{-yk) , 

where \arg{—yk) \ = | , is helpful. 
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Appendix H. Scattering solutions 

The usual scattering solutions of Lxijj{x,p) = for continuous real energies E = 
p^/ 2m > may be characterized according to the homogeneous BC they obey. 
Contrary to the resonant solutions, which are subject to two BCs, each one of the 
scattering solutions solution is subject to one homogeneous BC, alternatively at x = 
or at X = L : 

-ipC(o) 

-ip^riL), (H.i) 



5x.0r*(^)io = 



where r{l) labels the right (left)-moving character of the impinging wave. Notice that the 
reverse notation is often found in the literature, with r{l) indicating the wave coming 
from the right (left) of the barrier. 
These solutions, for instance 

= [ 0{-x){e'P^+Re-'P'')+9{x)9{L-x){P e'P'^'+Q e-'P''')+9{x-L) T e'P''] , 

are (5-normalized in energy. 
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